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ABSTRACT 

■ We suggest that some observational features of high-energy radiation from pulsars should be explained in terms of 
three dimensional geometric models, e.g. the phase-resolved X-ray and 7-ray spectra and the energy dependent light 
curves from various pulsars. In this paper, we present a three dimensional pulsar outer-magnetospheric gap model to 
explain these observational features. The outer-magnetospheric gap is proposed to form near the null charged surface 
and extend toward the light cylinder. The other two geometric dimensions of the outer-magnetospheric gap, i.e. the 
vertical size and the azimuthal extension can be determined self-consistently. We apply this model to explain the 
observed phase-dependent spectra and the energy-dependent light curves of various pulsars. 

INTRODUCTION 

The non-thermal high-energy radiation from rotation-powered pulsars is believed to be emitted from their mag- 
netosphere. In the past one decade, there has been a lot of progress in detecting and understanding high-energy 
radiation from rotation powered pulsars. In particular, the X-ray and 7-ray data obtained by the satellite observatories 
ROSAT, ASCA, RXTE, BeppoSAX, RXTE, CGRO, Chandra and XMM-Newton provide very important constraints 
for the theoretical models, which are used to explain the high-energy radiation from rotation powered pulsars. The 
observed data is so rich that the local properties in the magnetosphere, e.g. retarded relativistic effect, the time of 

■ flight, the strength of the local magnetic field, the local radius of curvature, the local soft photon density etc., are nec- 
i essarily to be used in order to explain the observed phase-dependent spectra and energy-dependent light curves from 

pulsars. These local properties associated with the observed photons not only depend on the global pulsar parameters, 
i.e. rotation period and surface magnetic field, but also depend on the angle between the magnetic axis and the rotation 
axis (the inclination angle a) as well as the viewing angle {Q. 

EGRET has observed the phase-resolved emission characteristics such as pulse profiles and phase-resolved spec- 
tra of Crab, Vela and Geminga Pulsars (Thompson et al. 1996). Fierro et al. (1998) have divided the observed photons 
of the Crab pulsar in the energy range from 100 MeV to 10 GeV into eight different phase bins and have shown that 
the spectra of each phase bin is very different. If a power-law fitting is used, the spectral indices of these phase bin 
can vary from -1.71 to -2.65. The observations of X-ray emission from the Crab pulsar also indicate that its spectral 
energy distribution is changing within its double peak pulse profile. Pravdo et al. (1997) analyzed the RXTE data 
between 5 and 250 keV, and showed an evolution of the spectral index across the X-ray pulse in a reverse S shape. 
Massaro et al. (2000) presented the phase-resolved analysis results based on BeppoSAX data. Kuiper et al. (2001) 
have summarized the basic observational properties of the Crab pulsar from soft X-rays to high-energy gamma-rays. 
Here, we summarize the basic observed properties at X-rays are as follows: (i) a double pulse profile with a bridge 
separated by ~ 0.4 in phase has been observed; (ii) the spectra of the two peaks soften and the softest spectrum is 
the first peak; (iii) the spectrum the peaks hardens; and (iv) spectral indices are clearly increasing with energy over 
all the phase interval. The Chandra X-Ray Observatory has also observed the Crab pulsar using the Low-Energy 
Transmission Grating with the High-Resolution Camera. Time-resolved zeroth-order images reveal the pulsar emits 
X-rays at all pulse phases (Tennant et al. 2001). A preliminary analysis of the dispersed data indicates that the spectral 



indices evolve as a function of pulse phase (Weisskopf, 2002). The Chandra result is in semi-quantitative agreement 
with previous measurements (e.g. Pravdo et al. 1997; Massaro et al. 2000) at various energies. However, the Chan- 
dra results extend the phase coverage through pulse minimum. In summary, observations strongly suggest that these 
phase-dependent properties come from different parts of the pulsar magnetosphere. The question is how to explain 
theoretically the phase dependent properties of the observed high-energy radiation from the Crab pulsar as well as 
from other gamma-ray pulsars. We suggest that the above examples indicate the necessity of using three dimensional 
outer magnetospheric gap (hereafter outer gap) models to explain the observed data. 

THREE DIMENSIONAL OUTER GAP MODELS 

The outer gaps, powerful acceleration regions, can form in the vicinity of "null charge surface" (17 • B = 0) 
(HoUoway, 1973; Cheng et al. 1976) because the charged carriers on each side of the null charge surface have opposite 
charges. In fact, the charge density of the magnetosphere in the corotating frame of a neutron star is (Goldreich and 
Julian, 1969) po '-^ —{ft ■ B/27rc), where B and ft are the magnetic field and angular velocity of the neutron star. 
The charge density will change sign when a global current flows through the null surface where O • B = 0. As 
a result, a charge-deficient region (p « 0) in the outer magnetosphere near the null surface will be formed. Any 
deviation of the charge density from po results in an electric field along B. Cheng et al. (1986a,b, hereafter CHRI and 
CHRn respectively) argue that this electric field can become strong enough to accelerate pairs to ultra relativistic 
energies. These pairs could radiate 7-ray tangential to the curved B field lines. These "curvature 7-rays" are 
further converted into pairs via 7 + 7 e"*" + e^. Therefore, in order to keep a steady state current flow and 
the charge density po in the regions outside the gap, the gap will grow until it is large enough and the electric field is 
strong enough to maintain a copious supply of charges to the rest of the open field line region. If the gap ends in a 
region pQ 7^ 0, charges from the surrounding region will flow in through the end. If both ends are located on the null 
surface, any produced in the gap will act to replace the charge deficiency inside the gap, and finally the gap will 
be filled up. However, if a vacuum gap extends to the fight cyfinder, charged particles created in the gaps will escape 
from the magnetosphere, so the gap will not be quenched. Hence, stable outer gaps (if they exist) are those from the 
null surface to the light cylinder along the last closed field lines. In each outer gap, the inner boundary of the outer 
gap lies near the intersection of the null surface where p = and the boundary of the closed field lines of the star on 
which the magnetosphere current does not flow. The thickness of the outer gap is bounded from above by a layer of 
electric current which contributes a surface charge density. 

According to CHR model, four outer gaps exist in the open zone in the plane of (ft, p) (two of them are 
topologically connected in three dimensional space), but only two longer outer gaps should give observed fan beams. 
They argued that these two longer outer gaps may create enough 7-ray and pairs to quench the two shorter, less 
powerful ones. In the CHR model, the 7-ray emission is approximated to occur only along the last closed field fine 
in the plane of the dipole and rotation axes. Because the charged particles of both positive and negative charges are 
accelerated in the gap, the emissions should beam both toward and away from the pulsar. Therefore, the observed 
fan beams consist of those coming from different gaps and the measured phase separation between the two peaks 
is determined by the time travel difference between these two outer gaps, relativistic aberration of emission and the 
bending of the magnetic field lines near the light cylinder. The emission from each peak is highly cusped because of 
the relativistic aberration, so there will be some bridge emission but very little other offpulse emission. Obviously, 
the pulse profiles of CHR model are not consistent with the observed those. So the three dimensional description of 
the outer gaps is necessary. 

After studying the 7-ray production and light curves for various magnetosphere geometries based on the CHR 
model, Chiang and Romani (1992) assumed that gap-type regions could be supported along all field lines which 
define the boundary between the closed region and open field line region rather than just on the bundle of field lines 
lying in the plane containing the rotation and magnetic dipole axes. In this case, photons are generated which travel 
tangential to the local magnetic field lines, and there are beams in both the outward (away from the neutron star) 
and inward directions, because the accelerating gaps are populated by pair production. They considered the pulse 
profile of radiation produced in the outer gap and showed that a single pole will produce a broad, irregularly-shaped, 
emission which is particularly dense near the edge. As a result, double 7-ray pulses will be observed when the fine of 
sight from the Earth crosses these enhanced regions of the 7-ray beam, while the inner region of the beam provides a 
significant amount of emission between the pulses. With a proper choice of the observer viewing angle, a wide range 
of peak phase separations can be accommodated. Furthermore, Chiang and Romani (1994) refined the calculation of 



high-energy emission from the rotation-powered pulsars based on the CHR model. Their major refinements include 
(i) the approximate location of the emission at each point in phase along a given line of sight was inferred by using 
a pulse phase map, and (ii) because the spectral emissivities at different emission points will differ, so the outer 
gap is divided into small subzones in the plane containing the rotation and dipole axes. The photon densities and 
beaming directions for different zone are also different, in which case the particle transport needs to be considered. 
Under their refinements, they found that the spectral variation of the 7-radiation over the pulsar period is the result 
of the different emission processes which play a role throughout the outer magnetosphere, however, they were not 
able to obtain a self-consistent spectrum which resembled the observed high-energy spectra, and they attributed this 
shortcoming to the inability to model appropriately the extremely complex emission processes and their interactions. 
Subsequently, Romani and his co-workers (Romani and Yadigaroglu, 1995; Yadigaroglu and Romani, 1995; Romani, 
1996) have improved their three dimensional models and successfully explain the high-energy emission features of 
pulsars including the phase-resolved spectra of the Vela pulsar. However, in their model they have assumed that there 
is only one single outer gap and only outgoing current. These two assumptions do not have real physical justification. 

CRZ MODEL 

Cheng et al. (2000 hereafter CRZ) re-consider the three dimensional magnetospheric outer gap model, following 
the important ground-breaking work of Romani and co-workers. But instead of assuming a single outer gap with 
only an outgoing current, and no restriction on azimuthal directions, they use various physical processes (including 
pair production which depends sensitively on the local electric field and the local radius of curvature, surface field 
structure, reflection of pairs because of mirroring and resonant scattering) to determine the three-dimensional 
geometry of the outer gap. In their model, two outer gaps and both outgoing and incoming currents are in principle 
allowed, but it turns out that outgoing currents dominate the emitted radiation intensities. Furthermore, the three 
dimensional structure of outer gaps is completely determined by pair production conditions. Since the potential drop 
of the gap is AF 6.6 x 10^^ f^BuP''^ Volts, where P = 27r/J7 is the rotation period, CI is the rotation angular 
velocity, B12 is the surface magnetic field in units of 10^^ Gauss, /o = h{< r >)/Rl, h{< r >) is the average 
vertical separation of the gap boundaries in the (CI, ji) plane and Rl = c/fi is the light cylinder radius, and < r > 
is the average distance to the gap; its value depends on magnetic inclination angle a{< r >~ Rl/'^)- The particle 
current passing through the gap is Ngap = 3 x 10^°/o^Si2P~^ s~^, where ^ = A$/27r; A$ is the transverse 
((?!)-direction) extension of the gap. Each of the charged particles inside the gap will radiate high-energy curvature 
photons with a characteristic energy E^^ (/o) = 2 x ^^i?^2^P' ''/^ eV. About half of Ngap will move toward the 
star. Although they continue to radiate their energies on the way to the star, they still carry 10.5P^/^ ergs of energy 
on to the stellar surface. The energy will be radiated back out in hard X-rays. However, resonant scattering with pairs 
near the star may reflect hard X-rays back to the stellar surface (Cheng et al. 1998; Wang et al. 1998), to be re-emitted 
as soft X-rays with a temperature ^ 3.8 x IQ^ fl^'^^^l'^Bl!^^ p-^l^"^ K. 

The X-ray photon density is very low but each pair produced by an X-ray-curvature photon collision in the outer 
gap will emit almost 10^ curvature 7-rays for further pair creation in that gap. Once the pair production threshold 
condition kTgE^ > (mgC^)^ is satisfied, the gap is unlikely to grow much larger. This pair production condition 

gives /o = 5.5P26/2l5-4/7^1/7_ 

Here, ^ is still an unknown quantity. However, /o is weakly dependent on ^ which 

is likely of order of unity. In first approximation, they assume /o = h.bP^^/^^ Bi^^^ (Zhang and Cheng, 1997). 
To determine ^, they consider local pair production processes. The pair production per unit length inside the gap 
is a decreasing function of r. According to CHR model, oc r~^/^ for the thin outer gap (e.g. the Crab pulsar), 
which gives E^{r) oc r~^/^ after using the large r hmit s{r) = {rRi)^^"^. Since E^ is only weakly dependent on r, 
they assume cr^-y ^ const. The local pair production per unit length is N^±{r) = (1 — e~'^'i"')N^{r) « Tj^Nj{r), 
where Tjy = nx{r)(7yyl{r) is the local optical depth, nx = R^Tfa /r'^kTsC is the X-ray number density at r, 
l{r) (2s(r)/(r)i?L)^/^ is the local optical path, /(r) = h{r)/RL is the local vertical extension of the gap (since 
B{r)h?{r) is a constant, which gives /(r) oc r^/^ and /o ~ /{Rl/^) ), and A^-y = eE^^{r)/E^{r) is the number of 
curvature photons emitted at r per e+/e~ per unit length. Then N^±{r) oc r~^^/^. Since most pairs are produced 
near the null surface where r = rj„, so the pair production take place mainly in the range rin < r < rum where 
ru„, is estimated as ru.„,N^±{rii,n)/rinN^±{rin) ~ {runi/rin)~^/^ ~ 1/2, which gives r;i„, 6ri„. This hmits pair 
production both along the field lines and in transverse directions, and gives A<i? '-^ 160° by using the parameters of 
the Crab pulsar. 



Within the pair production regions, outgoing and incoming directions for particle flows are allowed. For r > 
only outgoing current is possible. Figure 1 shows our 3D outer gap structure. Two pencil beams represent radio beams. 
The light budge grey shadow is the surface of last closed magnetic field fines, two dark surfaces represent the upper 
boundaries of two outer gaps. So the structure of the outer gap starts from the null surface and end at the light cylinder. 
The lower boundary is last closed field surface and the upper boudary is shown in Figure 1. 




Fig. 1 . 3D structure of the outer-magnetospheric 
gaps. 



Fig. 2. Emission projection onto tlie (C, plane 
and pulse profile. The emission consists of the 
emission outwards from both outer gaps and in- 
wards only from the region {rum-rin) of both outer 
gaps. The outer gaps are limited along the az- 
imuthal direction by pair production. Crab param- 
eters, a = 65°, C = 82° and A$ = 160° are used. 



EMISSION MORPHOLOGIES AND LIGHT CURVES 

In this section, we discuss the morphological features of emission from the three dimensional outer gaps. It 
is important to note that in order to properly describe the three-dimensional pulsar magnetosphere three-dimensional 
rotating dipolar magnetic field must be used. Chiang and Romani (1994) and Romani and Yadigaroglu (1995) assumed 
the entire locus of points on the last closed surface bounded by the null surface and the light cylinder as giving emission 
in an outer gap. In our model, the emission-producing outer gaps are limited along both radial and the (/)-directions. 
For example, the extension of the outer gap on the ^-direction is about 160° for Crab parameters with a = 65°. We 
describe the new photon emission morphologies below. 

High-energy photons will be emitted nearly tangent to the magnetic field lines in the corotating frame because 
of the relativistic I/7 beaming inherent in high energy processes unless |E x B| ~ B^. Then following Romani and 
Yadigaroglu (1995), we assume relativistic charged particles in the open zone radiate in their direction of propagation, 
i.e. along the magnetic field lines in the corotating frame. For each location within the open zone the direction 
of emission expressed as (C,^) is calculated, where ^ is the polar angle from the rotation axis and $ is the phase 



of rotation of the star. Effects of the time of flight and aberration are taken into account. A photon with velocity 
u = {ux, Uy,Uz) along a magnetic field line with a relativistic addition of velocity along the azimuthal angle gives an 
aberrated emission direction u' = (u^, Uy,u'^). The time of flight gives a change of the phase of the rotation of the 
star. Combining these two effects, and choosing $ = for radiation in the (x,z) plane from the center of the star, 
and $ are given by cosC = u'^ and $ = —(/)u' — f • u' , where (pu' is the azimuthal angle of u' and r is the emitting 
location in units of Rl- We project photon emissions on the (C, plane and observe the emission patterns on the 
sky. In {(,, $) plane, the null surface can be determined easily because it consists of the points at which magnetic 
field Unes are perpendicular to the rotation axis. For a field line, the null charge crossing is where the projected line 
crosses the equatorial line {C, = 90°). As an example, we show the projections of photon emission both outwards and 
inwards from one outer gap in {C„ ^) plane for the magnetic inclination angle (q = 65°) in the upper panel of Figure 
2, where the extension of the outer gaps on the ^-direction is assumed to be 180° and {rum — Tin) ~ 0.55i?i:(CRZ). 
As mentioned above, the pulse profile depends on emission location and viewing angle. Here we use Crab parameters 
and the viewing angle C, = 82° ( the dashed line in panel A of Figure 2) to construct the light curve. The intensity of 
the light curve at each phase (panel B of Figure 2) is proportional to the number of interceptions between the dashed 
line and the sohd lines in panel A of Figure 2. 

PHASE-RESOLVED SPECTRA OF THE CRAB PULSAR 

In this section, we describe how to calculate the phase-dependent spectra of pulsars. We will use the parameters 
of the Crab pulsar as example. Because the Crab pulsar outer gaps are thin, we use the electric field of the CHR 
model: 

nB{r)a\r) mr)f{r)Rl 

-^iiv) = r\ — ^ r\ ' 

" cs[r) cs(r) 

where a(r) is the thickness of the outer gap at position r, the radius of the curvature s (ri?^)^/^, and /(r) = 
a{r)/RL is the local fractional size of the outer gap. Since the magnetic flux subtended in the outer gap should be 
constant in the assumed steady state, 

/ „ x3/2 

/w~ (^) . (2) 

where /{Rl) is estimated by the pair production condition described in previous section(N.B. There are two possible 
ways to determine /{Rl), Eq. (6.7) of CHR II or Eq. (22) of Zhang and Cheng (1997), but they come out very close 
to each other.). The local Lorentz factor of the accelerated electrons/positrons in the outer gap is 

^e{r) = i^-^^^eE^\{r)cj . (3) 

Because of the high soft photon density, the high-energy emission from the Crab pulsar is described by synchrotron 
self-Compton process. 

In an outer gap, the number of primary charged particles in a volume element AV is roughly given by dN = 
nojAAAl, where noj = ^ ■ B/27rec is the local Goldreich-JuUan number density, BAA is the magnetic flux 
through the accelerator and Al is the path length along its magnetic field Unes. Using the thin gap approximation, the 
total number of charged particles in the outer gap is 

N ~ Rl, (4) 

47rce 

where $ ~ f {Rl)B{Rl)R\A(I) and Acf) is the angular range of the outer gap extending along the azimuthal direction, 
estimated in previous section. These primary pairs will lose their energy by radiating curvature photons with a 
characteristic energy 

Ecurir) = \Hir)^y (5) 
2 s[r) 

The power into curvature radiation for dN pairs through in A unit volume is 



(6) 



where Icur = (r)c, is the local power into the curvature radiation from a single electron/positron. The spectrum 
of primary photons from a unit volume is 

d^N IcurriGJ 1 



(7) 



where Ey < E^ur- These primary curvature photons collide with the soft photons produced by synchrotron radiation 
of the secondary pairs, and produce the secondary pairs by photon-photon pair production. Although pair 
production inside an outer gap is limited to a small region (ri„ < r < rum), pair production outside the outer gap 
can cover a much wider range because the synchrotron photons produced by the secondary pairs are more abundant 
than the thermal photons from the stellar surface. The former cannot get into the outer gap because of the field line 
curvature (cf. CHR I) but they can convert most curvature photons from the outer gap into the secondary pairs. In a 
steady state the distribution of secondary electrons/positrons in a unit volume 

(fN 1 f d^N{E'y = 2E',) llcurncj, f E,ur\ 

dEe^ -z 7=, In , (8) 



dVdE^ J dVdEy £^c«r V 

with Ef. the electron energy loss into synchrotron radiation, 

2efB2(r)siii^3(r) ^ ^, ^ ^ 
3 rr?(? 

B(r) is the local magnetic field and /?(r) the local pitch angle, 



^e = -o ^23 Ml ' (9) 



sm/3{r) ^ sin /3{RL)i^—j . (10) 

sin (3{Rl) is the pitch angle at the hght cylinder. Then the energy distribution of the secondary electrons/positrons in 
volume AV{r), 

(d_N^^ ^ _^^y^^^ ^ ^ k^rngjAVir) ( E^^ 

V dEe J dVdEe ^ ' E, Ecur \ J 



The corresponding photon spectrum of the synchrotron radiation is 

!l/2,-3 

J->\ i ) ^^^^ fJ 

\ Je, 

where x = Ej/Egyn 



._ 3V^:^I?(r) sin/3 1 f""-^ (dNir)\ 



3 / EeY heB{r) sin I3{r) 
2 \ mc^ J mc 



is the typical photon energy, and F{x) = x Kv:,i^{y)dy, where K^/^{y) is the modified Bessel function of order 
5/3. Similarly, the spectrum of inverse Compton scattered photons in the volume AV{r) is 

d^N{r)ics _ n ^^^^^^^ ^^p^^^ p^^^^^ 



where 



dE^dt J 



and 



3aTC 

F[e, Ey, Eg 



A{Ee/mc^Y e 



2glng + (l + 2g)(l-g)+ ^^^^|'|^^J^ 



(16) 



with r = 4e{Ee/mc^) /mc^, q = Ei/T{l - Ei) with Ei = Ey/Ee and l/4:{Ee/mc^) <q<l. The number density 
of the synchrotron photons with energy e is 



where F^yn is the calculated synchrotron radiation flux, and Ari is the usual beam solid angle. 

Since the outer gap of the Crab pulsar is very thin, it is sufficient to use one representative surface to calculate 
the high-energy radiation. For a given viewing angle not only the Ught curve can be determined but also the exact 

emission regions in the outer gap are known. Figure 3 shows the emission trajectories in the outer gap. Once the radial 
distances of the emission regions are determined, the spectrum of photon emission can be calculated for a given radial 
distance r. Figure 4 shows the comparison between the observed phase-resolved spectra of the Crab pulsar and the 
calculated spectra in the energy range of 10 MeV to 10 GeV (Zhang et al. 2000). The more detailed phase-resolved 
spectra only in the energy range of EGRET can be found in CRZ. Figure 5 show the model energy dependent light 
curves of X-rays in four different energy channels (Zhang and Cheng, 2001b). Figure 6 shows the phase dependent 
spectral indices of X-rays in four different energy ranges. 



PI TW1 Bridge LW2 P2 




50 70 90 110 130 150 170 190 

Phase (Deg.) 

Fig. 3. The variation of radial distance witli pulse phase for the Crab. The inclination angle is 65°. Five 
regions for different pulse phase are indicated. 



THE PHASE-RESOLVED SPECTRA OF THE GEMINGA PULSAR 

Fierro et al. (1998) have shown the observed light curve and phase-resolved spectra of high-energy 7-rays of 
Geminga pulsar detected by EGRET. The observed pulse profile by EGRET indicates that the phase separation is 
0.49 lb 0.05. In order to obtain the observed phase-resolved spectra, Fierro et al. (1998) divided Geminga pulsar phase 
into 8 parts: leading wing 1 (LWl), peak 1 (PI), trailing wing 1 (TWl), bridge, leading wing 2 (LW2), peak 2 (P2), 
trailing wing 2 (TW2) and offpulse (OP). The phase intervals widths of all these parts are 0.11, 0.09, 0.11, 0.15, 0.13, 
0.13, 0.08 and 0.21 respectively. They have obtained the spectra for these different phase intervals and shown that the 
spectral indices change as a function of phase. 

Now we apply to this model to explain the pulse profile and phase-resolved spectra of the Geminga pulsar. The 
key differences between the Geminga pulsar and the Crab pulsar are: (1) According to the model described in previous 
section, we can find that the fractional height of the outer gap for Geminga is /o 0.7, which is a very thick gap. 
Certainly we cannot use a single surface to represent the emission regions. We have used five layers to approximate 
the emission regions (cf. Figure 7a). (2) The mean free path of the primary photons from the gap is longer than the 
light cylinder. So the observed 7-rays are curvature photons from primary charged particles instead of synchrotron 
photons from the secondary pairs like in the case of the Crab pulsar. In order to calculate the light curve and phase- 
resolved spectra, we also need to know a and Since the Geminga pulsar is a radio-quiet pulsar, these two parameters 
are difficult to know. Cheng and Zhang (1999) proposed a model for X-ray emission from rotation-powered pulsars. 
They applied this model to the Geminga pulsar and found that the magnetic inclination angle is ~ 50°. Here, we use 
this value of the magnetic inclination angle. Furthermore the 7-ray light curve of Geminga is nearly ~ 180° so the 
viewing angle must be very close to ~ 90°, which is consistent to the fact that its radio beam cannot be within the line 
of sight. In Figure 7b, we have compared the model phase-resolved spectra and the observed data. The more detailed 
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Fig. 4. Phase-resolved 7-ray spectra from lOMeV to lOGeV for peak 1 , bridge, peak 2 and pliase-average 
of the Crab pulsar. Observed data are taken from Ulmer et al. (1995). 



calculations on the X-ray and 7-ray light curves of the geminga pulsar can be found in Zhang and Cheng (2001a). 

APPLICATIONS TO OTHER CRAB-LIKE PULSARS: PSR B0540-69 AND PSR B1509-58 

We have also applied this three dimensional outer gap model to other Crab-like pulsars:PSR B1509-58 and PSR 
B0540-69, and model the light curves and the spectra of optical, X-rays and 7-rays from these two pulsars (Zhang and 
Cheng, 2001b). Although the mean free path of the primary photons from the outer gaps of these two pulsars is much 
shorter than the light cylinder and hence the radiation from secondary pairs still dominate in the observed spectrum, 
the thickness of the outer gap in these two pulsars are 0.25 and 0.3 respectively. Again the emission regions cannot be 
approximated by a single surface. In Figure 8, we have compared the model spectra and the observed spectra of PSR 
B0540-69 from optical to gamma-rays. The more detailed comparison between model results and the observed data 
of PSR B 1509-58 and PSR B0540-69 can be found in Zhang and Cheng (2001b) 
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Fig. 5. Expected X-ray pulse profiles of Crab pul- 
sar for four different energy bands. The first model 

light curve in the energy band from 0.3 - 4.2 keV Fig. 6. Comparison of expected spectral indices 
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magnetic inclination and viewing angles are as- from Massaro et al. (2000) and Weisskopf (2002). 
sumed to be 65° and 82°. 



SUMMARY AND DISCUSSION 

We use a 3-D model magnetosphere to model the observed light curve and the phase-resolved spectra of the 
7-ray pulsars. In our model, the local photon-photon pair production in the outer gaps limits the extension of the outer 
gaps along the azimuthal direction. We find that the two topological disconnected outer gaps, with some extension 
along the azimuthal {(/)) direction, exist in the pulsar magnetosphere. Double-peaked pulse profiles with varying phase 
separation, depending on viewing angle, and strong bridge emission occur naturally, as in the single pole outer gap 
model. In case of the Crab pulsar, pair production is not limited to inside the outer gap because the intense X-rays 
produced by secondary pairs in the outer-magnetosphere of the Crab pulsar can convert most curvature photons into 
pairs outside the gap. The observed spectrum of the Crab pulsar results from a synchrotron-self-Compton mechanism. 
We (Cheng et al. 2000) have tried to apply the 3D outer gap model to explain the Crab pulsar's phase-resolved spectra 
and find some discrepancies between model results and the observed data. One possibihty is that the data presented 
by Fierro et al. (1998) contains the DC emission from the Crab/Nebula. 

Although our three dimensional model can explain the observed phase-resolved spectra and the energy dependent 
light curves of various gamma-ray pulsars including the Crab pulsar and the Geminga pulsar reasonably well, the 
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Fig. 7a. The variation of radial distance with the pulse phase for different outer gap surfaces of the Geminga 
pulsar. The inclination angle is 50°. Five regions for different pulse phases which are the same as observed 
those are indicated. 



present model is subjected to number of limitations. First, the photons in the trailing wing 1, off -pulse and leading 
wing 1 in the light curve of the Crab pulsar (Fierro et al. 1998) cannot be explained by this model. Secondly, although 
the model Chandra light curve of the Crab pulsar and the observed one (cf. Figure 4 of Weisskopf, 2002) are very 
similar, we can only fit the phase-dependent spectral indices in half of the period (cf. Figure 6). Again, we can explain 
the spectra of PSR B0540-69 and PSR B 1509-58 reasonably well but the model X-ray light curves are narrower than 
the observed ones. 

In order to explain these discrepancies between the model results and the observed data, some improvements 
can be made in our model. (1) In our calculations, photons are assumed to be emitted tangent to the local field lines. 
Therefore, the emission light curves (cf. Figure 3) have sharp edges in both ends. In fact, some of secondary pairs 
have quite large pitch angles whose emission will not be tangent to the local field lines. Also for those pairs created 
near the null charged surface even they begins with small pitch angles, but when they stream towards the star and part 
of them will reflect back outwards due to the magnetic mirroring effect (Ho, 1988). Then they could end up with very 
large pitch angles. (2) The realistic magnetic field configuration in the pulsar magnetosphere could be different from 
a simply rotating dipolar magnetic field structure. (3) The particle energy density and the magnetic energy density 
are comparable near the light cylinder. Charged particles are not necessary straightly moving along the field lines. 
These pairs will have a much larger pitch angles and their radiation should not be restricted between pulses. (4) The 
charged current inside the outer-magnetospheric gap is assumed to be Goldreich- Julian current (1969), which must 
be the maximum value. In order to determine the real current flow inside the gap, it is necessary to solve a more 
consistent electro-dynamic model (Hirotani, 2002). 

Finally, we would like to make a few remarks. The better test of the three dimensional outer gap model is 
to compare the model results with the phase-resolved spectrum from X-ray band to 7-ray band together, instead of 
comparing them separately. For the Crab pulsar, Kuiper et al. (2001) have published a full observational (phase- 
resolved) picture from soft X-rays up to high-energy gamma-rays. For the Geminga Pulsar, Jackson et al. (2002) 
have presented the combined data of ASCA, CGRO and RXTE. For other Crab-like pulsars, e.g. PSR B 1509-58 ( 
Kuiper et al. 1999; Cusumano et al. 2001) and PSR B0540-69 (de Plaa et al. 2003), more X-ray and 7-ray have been 
reported. We must carefully analyze these new data to see if they really support a simple three dimensional outer gap 
model. In fact, de Plaa et al. (2003) have point out that there are noticeable discrepancy between latest data and model 
predictions (Cheng et al. 2000). They suggest that these discrepancies may result from the uncertainties in the pulsar 
geometry. 
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